Pair distribution function and structure factor of spherical particles 



O 
O 



CO 



Rafael C. HowelQ 

Materials Science and Technology Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 

Thomas ProfferQ 

Lujan Neutron Scattering Center, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 

Steven D. Conr adsorfj] 

Materials Science and Technology Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 

(Dated: February 2, 2008) 

The availability of neutron spallation-source instruments that provide total scattering powder 
diffraction has led to an increased application of real-space structure analysis using the pair distri- 
bution function. Currently, the analytical treatment of finite size effects within pair distribution 
refinement procedures is limited. To that end, an envelope function is derived which transforms the 
pair distribution function of an infinite solid into that of a spherical particle with the same crystal 
structure. Distributions of particle sizes are then considered, and the associated envelope function is 
used to predict the particle size distribution of an experimental sample of gold nanoparticles from its 
pair distribution function alone. Finally, complementing the wealth of existing diffraction analysis, 
the peak broadening for the structure factor of spherical particles, expressed as a convolution derived 
from the envelope functions, is calculated exactly for all particle size distributions considered, and 
peak maxima, offsets, and asymmetries are discussed. 
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I. INTRODUCTION 



There currently exists much scientific interest of the 
physical and chemical properties of nanoparticles and 
nanodomains. Generally investigations of particles with 
diameters of the order of one micron are limited to pow- 
ders, and from the diffraction data one performs Rietveld 
refinement to extrapolate their structure and size, by 
using the Debye-Scherrer formula^ extensions such as 
Debye function analysis^ and small-angle scattering. 3 
Alternatively one derives the pair distribution func- 
tion (PDF) from the diffraction data^ thus facilitating 
the study of size, correlated atomic motion^ short- to 
medium-range order ^ and other phenomena more ap- 
parent with a real-space treatment. 

A number of techniques have been developed to infer 
the size and structure of particles from the PDF alone. 
For example, the local structure has been determined by 
an interpretation of the first peaks of the PDF^i and the 
size has been estimated from a Fourier transform of the 
wide-angle Debye-Scherrer diffraction pattern^ In this 
manuscript, a rigorous approach to the determination of 
particle size is taken by rederiving the pair distribution 
function of a single spherical particle, expressed as an 
envelope function that multiplies the PDF of an infinite 
crystal with the same crystal structured A general class 
of distributions of particle sizes is then considered, and 
an associated distributed envelope function is obtained. 
Using experimental PDFs of both bulk gold and gold 
nanoparticles, calculated from high-Q ncutron-powder- 
diffraction data, a distributed envelope function is used 
to transform the PDF of bulk gold to give a best fit repli- 



cation of the PDF of gold nanoparticles. Based on the 
parameters of this envelope function, the particle size 
distribution of the gold nanoparticles is predicted and 
compared to that obtained experimentally. Finally, a re- 
lationship between the real-space envelope function and 
a Q-space convolution function is established, the latter 
of which is to be applied to the structure factor of an 
infinite crystal to obtain that of a distribution of spheri- 
cal particle sizes. The analytical form of the convolution 
function allows for a quantitative analysis of Bragg peak 
maxima, widths, and asymmetries as a function of peak 
position and particle size distribution. A thorough ac- 
count of the relationship between the PDF and structure 
factor can be found in the literature i&iSi 

The formalism used is identical to Peterson et al.^ 
with comparisons to other definitions and nomenclature 
found in KeeniiS With little effort, the conclusions made 
here regarding the PDF and structure factor of spherical 
particles can also be carried over to embedded spheri- 
cal domains, with the stipulation that the individual do- 
mains be uncorrelated with each other (i.e., they have 
random orientations), and uncorrelated with the host 
matrix. Indeed, diffraction analysis using spherical ge- 
ometries has already been succcessful in studies of water 
in mesopores and micropores^ 3 - as well as distributions 
of particle and void sizes in NMR cryoporometryii^ii&i 
The structure factor analysis of spherical domains with 
well-defined atomic structures fits within the more gen- 
eral context of disorder within crystals, thus contribut- 
ing to the analytical treatment of the associated diffuse 



scattering 
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It also provides a means of quantifying 



the diffraction limit with respect to localized lattice dis- 
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tortions, with spherical domains being a special case to 
be considered within the large class of nanoscale hetero- 
geneities already studied^ 



II. THE PAIR DISTRIBUTION FUNCTION OF 
A SINGLE SPHERICAL PARTICLE 



The microscopic pair density gives a distribution of 
atomic pair distances r in a sample, weighted by the 
pair's scattering lengths 



p(r) 
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The following method of constructing p(r) will be useful 
in determining its form with respect to spherical parti- 
cles. For each r, define a spherical shell with this radius. 
Let the center of this sphere coincide with the position of 
an atom i, and record as weighted <5 functions every atom 
j that intersects the spherical shell. Finally, divide the 
result by its surface area and the total number of atoms 
N. The calculation of p(r) for a spherical particle lim- 
its the position of the center of the sphere to the atoms 
within the particle itself, whereas the sphere's surface can 
extend beyond. Note that p(r) for a single particle may 
only be a part of the total microscopic pair density of a 
solid or solution. 

The contribution to p(r) of all atomic pairs {ij} within 
a spherical particle is limited to the range < < 2R, 
where R is its radius. If the particle is in solution, or 
we consider an ensemble of identical particles with ran- 
dom orientations embedded within a host lattice, then 
p(r) = po for r > 2R, where po is the constant atomic 
number density outside the particle. To simplify the fol- 
lowing, let po be equal to the number density of the par- 
ticle itself, or take po = for empty space. The essence 
of the problem addressed hereafter is to quantify the rel- 
ative population of atomic pair distances between any 
two atoms within the particle and pair distances where 
one atom resides outside of the particle. That is, p(r) 
will have r-dependent contributions from both the mi- 
croscopic pair density p c (r) of an infinite crystal and the 
uncorrelated outside structure po, 



p(r)=f e (r,R) Pc (r) + [l-f e (r,R)]p , 



(2) 



where < f e (r,R) < 1. When R -> oo, f e {r,R) = 1 
for all r and p(r) is that of an infinite crystal. When 
R -> 0, f e (r,R) = for all r and p(r) = p . f e (r,R) 
is the envelope function that we now derive when R is 
between these limits. 

Consider a point within the particle whose position is 
given by the vector r', with the center of the particle 
defining the origin. Orient the point and the particle so 
that r' aligns with the z axis, as shown in Fig. ^ A 
spherical shell of radius r around this point will either 
be enclosed by the particle if r < R — r' , intersect the 




FIG. 1: A spherical particle with radius R. An atom a dis- 
tance r' from the center of the particle can have a shell of 
radius r that is only partially embedded within the particle. 



surface of the particle if R — r' < r < R + r' , or enclose 
the particle if r > R + r' . 

When R — r' < r < R + r' , a line from any point on the 
circle of intersection and the position r' will meet the z 
axis at an angle a. The fraction of the surface of radius 
r around this point that is enclosed within the particle is 

r-2-n 

f{r',r,R) = — L — I rd<t> I rsm9d9 
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cos a). 
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Using the two right triangles from Fig. ^ the angle a 
can be expressed as 



R 2 - r' 2 - 



cos a — 



2r'r 



(4) 



The contribution of all such spheres enclosed within the 
particle is obtained by integrating /(/, r, R) over the re- 
maining positions r' in the region occupied by the par- 
ticle, taking care to consider when the shell of radius r 
extends outside the region. Using Eq. for r < 2R, 
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Finally, dividing by the total particle volume gives the 
envelope function 



f e (r,d) = 



1 - 



3 r 1 /ry 
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Q(d-r), 
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where d = 2R is the particle diameter and <d(x) — 0(1) 
for negative (positive) x is the Heaviside step function. 
f e (r,d) and its derivative are continuous for all positive 
r. 

The PDF of the particle is related to the microscopic 
pair density 

G{r,d) = Airr [p(r) - p ] 

= f e {r,d)G c (r), (7) 

where G c (r) = 4irr [p c (r) — po) is the PDF of an infinite 
crystal with the same crystal structure as the particle. 

The fraction of atom pairs residing within a spherical 
region is obviously not a continuous function of the re- 
gion's size, considering the discrete nature of a crystal. 
Therefore, one expects the derivation above to become 
less accurate for smaller particle sizes. To test this, the 
exact calculation of G(r) [using Eq. ^ for spherical par- 
ticles with an ideal fee structure with lattice constant a 
was compared to G(r) obtained by applying Eq. [7| to 
the PDF of an infinite ideal fee structure. A particle was 
constructed with n shells of atoms around a central atom, 
giving it a diameter of d — \phia (some shells contain 
no atoms). While the difference in the nearest-neighbor 
peak heights between the two was 7% for d = 2a, the dif- 
ference between all peaks rapidly diminished to less than 
1% for d — 4^2 a, and became negligible thereafter. 

III. DISTRIBUTIONS OF SPHERICAL 
PARTICLE SIZES 

The envelope function f e (r, d) depends on a single- 
particle diameter d. From a distribution P(d') that de- 
fines an ensemble of particle diameters d! , a distributed 
envelope function /de(0 can be constructed by weight- 



As with the single particle distribution G(r, d, n) — 
/de {r, d, n)Gc(r). 

Figure |21 shows the single particle size envelope func- 
tion f e (r, d) (Eq. and the distributed envelope func- 
tions /de (r, d, n) for n — 3 and n — 9 (expressed using 
the average diameters d = AD and d = 10D, respec- 
tively). The associated distributions are shown in the 
inset, and the case of n = 100 is included to illustrate 
the trend of these distributions (the single particle size is 
a S function distribution, and is not shown). The width of 
the distributions clearly dictates the shape of /de( ? "i d, n). 
A broader distribution gives the resulting envelope func- 
tion a longer tail at large r at the expense of a steeper 



ing individual single particle envelope functions with this 
distribution, 

POO 

/deM= / f e (r,d')P(d')dd'. (8) 
Jo 

Consider the following normalized distribution: 

where n is a positive integer and D is a positive real num- 
ber. The average particle diameter d for this distribution 
is related to the parameters n and D, 

POO 

d = (d') = / d'P(d') dd' = (n+ 1)D, (10) 
Jo 

and the characteristic width of the distribution is 

a=y/(d^-<P = -^=. (11) 
Vn + 1 

Hereafter the average diameter d will be used to identify 
the associated distributed envelope functions and their 
behaviors, thus facilitating comparisons between those 
obtained from various distributions and the single parti- 
cle size envelope function. The parameter D is still used 
within function expressions, however, to maintain their 
simplicity. The two are always related by Eq. Hoi 

When deriving a distributed envelope function from 
Eqs. [S]and [5] alone, one recognizes that a weighted sum 
of these distributions (each with unique n and D) yields 
a weighted sum of individual distributed envelope func- 
tions, so that the results hereafter can be easily adapted 
to a variety of distributions. For n > 3, a closed form 
expression of the resulting envelope function is given by 



(12) 

I 

decline for small r. Note that a universal condition for 
all spherical particle size distributions with average di- 
ameter d is 

f°° 3 
/(r = 0,d) = l and J f(r,d)dr=-d. (13) 



IV. COMPARISON WITH EXPERIMENT 

By comparing an experimentally obtained PDF of a 
collection of spherical-like particles with the theoretical 
expressions just derived, it should be possible to pre- 
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FIG. 2: The single particle size envelope function and two 
examples of distributed envelope functions, with n = 3 and 
n = 9. The associated distributions P(d! /d) are shown in the 
inset along with n = 100 (dotted). 



diet the particle size distribution used in the experiment. 
Neutron data of a 2 g batch of capped gold nanopar- 
ticles and a bulk gold fee powder reference were col- 
lected at T — 15 K on the neutron powder diffractome- 
ter (NPDF) at the Lujan Center at Los Alamos National 
Laboratory^ The PDF of both were obtained, and the 
bulk gold underwent full profile structural refinements, 
using PDFFIT, 21 to account for a host of influences such 
as correlated and uncorrelated atomic motions and in- 
strument resolution. Multiplying the refined PDF of the 
bulk data Gt,(r) with an envelope function derived from 
either a single particle size f e (r,d) or a distribution of 
particle sizes /de(?™, d, n) should give a good representa- 
tion of the ideal PDF for spherical nanoparticles of fee 
gold. A simplification is made by assuming that the crys- 
tal structure within the bulk gold and gold nanoparticles 
is identical. This simplification will be addressed when 
the results of the PDF comparisons are discussed. 

Using the PDF data of gold nanoparticles G np (r) the 
best fit diameter (and exponent for the distribution) is 
obtained by minimizing the root mean square deviation 
Sq between this PDF and the transformed bulk PDF. 
The deviation is given by 



51 = L- 1 



[G np (r)-f(r,d)G b (r)] 2 dr, (14) 
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where the lower bound is chosen to ignore spurious os- 
cillations in the experimental PDF, the upper bound is 
determined from the scattering resolution, and L is the 
difference between these two. The integrations are per- 
formed numerically with a lattice spacing of 0.01 A. 

Table [I] compares the parameters from the best fit par- 
ticle size distribution and the single particle distribution. 



Distribution 


n 


d (A) 


a (A) 


So (A 2 ) 


Single 




28.75 




1.243 


P(d',d,n) 


13 


29.80 


7.96 


1.238 


Experiment 




35.48 


13.12 





TABLE I: A comparison of the parameters from the best fit 
particle size distributions and the experimentally determined 
distribution. 



The predicted distribution of particle sizes, with n = 13 
and d n — 29.80 A, gives the smallest absolute deviation 
from the experimental PDF, with Sq = 1.238 A~ 2 . 
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FIG. 3: An offset comparison of the refined experimental PDF 
from the bulk gold, the predicted PDF from the single particle 
size and the distribution of particle sizes, and the experimen- 
tal PDF from the gold nanoparticles. 

Figure shows an offset comparison between the ex- 
perimental PDF of the gold nanoparticles, the best fit 
PDFs from the single particle size and the distribution 
of particle sizes, and the refined PDF of the bulk gold 
used to obtain the fits. The envelope function for the 
distribution of particle sizes allows for correlations over 
longer atomic separations without changing the behavior 
of short-range correlations, as is evident in the figure, and 
is thus a better fit to the experimental PDF than that 
of the single particle size distribution. The PDF derived 
from a single envelope function always truncates to zero 
for r > d. 

As a further test of the theory of spherical particles 
developed so far, a comparison of the particle size distri- 
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FIG. 4: A comparison of the experimentally determined par- 
ticle size distribution (histogram), and the predicted particle 
size distribution from the best fit PDF. 



bution given by minimizing Eq. 1141 with the distribution 
obtained experimentally is now possible. Using a JEOL 
2010 TEM with point-to-point resolution of 1.9 A, the 
diameters of 148 gold nanoparticles within the experi- 
mental sample were measured, giving the distribution of 
particle sizes shown in Fig. ^(histogram). The average 
particle diameter from this distribution is d = 35.48 A, 
and the width of the distribution is a = 13.12 A. These 
values are given in TableQ] Also shown in the figure is the 
predicted theoretical particle size distribution (d = 29.80 
A, a = 7.96 A). Obviously the best fit PDF predicts a dis- 
tribution of particle sizes with a smaller average diameter 
and width. This can be attributed to the presence of pen- 
tagonal twinning seen within the gold nanoparticles (see 
Fig. 2 of Page et alm&). The characteristic length scale of 
the twinned regions defines the scale of long-range order 
(an fee structure in this case) , which affects the resulting 
PDF. Each object might be better described as a collec- 
tion of single domain crystallites instead of a spherical 
particle with a uniform structure. Eauation ll4l compares 
their PDF with the bulk gold PDF, in which twinning 
did not occur and the structure was uniform, and thus 
delivers a best fit envelope function that represents the 
size of the coherent fee structure. The counting method 
used to construct the histogram in Fig. 0] however, con- 
sidered only total particle size, which contributes to the 
discrepancy between the two distributions. 



V. THE STRUCTURE FACTOR OF 
SPHERICAL PARTICLES 

The total-scattering structure factor of spherical par- 
ticles S(Q, d) can be obtained directly from their PDF 



by a Sine transform 

poo 

Q [S{Q, d) - 1] = / G(r, d) sin(Qr) dr 
Jo 

f{r,d)G c (r)sm(Qr)dr. (15) 



Let S C (Q) be the structure factor of the associated in- 
finite crystal. f(r, d) is an envelope function from any 
distribution of particle sizes, and can be expressed as the 
inverse cosine transform of a function f(Q,d), 



f(r,d) 



f(Q, d) cos(Qr) dQ. 



(16) 



Since G c (r) is the inverse Sine transform of Q [S C (Q) — 1], 
Eq. 1151 can be written not only as a Sine transform of 
a product of two functions, but also as a convolution of 
their respective transforms 22 



'(S(Q,d)-l) = 



[f(\Q-Q'\,d)-f(Q + Q',d)) 



xQ'(5 c (Q')-l) dQ' 



1 f°° 

S(Q, d) = — J [f(Q - Q', d) - f{Q + Q', d)] 

xQ'S c (Q')dQ'. (17) 

The absolute value of the argument Q — Q 1 can be disre- 
garded if one recognizes the convolution function f(Q,d) 
as an even function of Q. 

For a single particle size distribution, the convolution 
function is 



f e (Q,d) = / f e {r,d) cos (Qr)dr 
h 

f n x {Qd) + 7t4? + ^)> (18) 



my 



{Qdf 



where n\(x) = —cos(x)/x 2 — sin(a;)/a; is a spherical 
Bessel function of the second kind. f e (Q,d) has a half- 
width at half- maximum of approximately Q = 3.48/d. 
This convolution function is very similar to the expres- 
sion used for the intensity of scattering from randomly 
oriented identical spherical particles, and is often used 
in small-angle scattering analysis^ 1 but differs from the 
sine function derived from a Fresnel construction of a 
spherical crystal that is often used to characterize pow- 
der diffraction lines. 

The convolution function for the distribution of parti- 
cle sizes defined by Eq. EJcan be obtained by noting that 
each term in Eq. 1121 has a cosine transform proportional 
to 
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k e~* cos(Qr) dr = — - V (-l)^ 2 ( k + ^ (QD) J '. (19) 
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FIG. 5: The convolution functions to be applied to the infinite 
crystal structure factor QS C (Q), for the single particle size and 
two distributions of sizes n — 3 and n — 9. 

Figure El shows the convolution functions for the single 
particle size and the two distributions of sizes n — 3 and 
n = 9, considered before (expressed using the average 
diameters d = AD and d — 10D, respectively). For the 
single particle size, shoulders appear after the primary 
peak, a feature attributed to the similarity between a 
periodic envelope function and a triangular wave. When 
distributions of particle sizes are considered, the shape is 
always Lorentzian-like, due to the exponential behavior 
of the associated distributed envelope function. 

A broader distribution gives the resulting convolution 
function a longer tail at large Q at the expense of a 
steeper decline for small Q. Note that a corollary to 
Eq. HSlis 



f(Q = 0,d) = -d and - 

8 7T 



f(Q,d)dQ = l. (20) 



To summarize, the structure factor of a distribution of 
particle sizes can obtained by convoluting the structure 
factor of the associated infinite crystal with an envelope 
convolution function f(Q, d). The result will be a broad- 
ening of the Bragg peaks of the infinite crystal, but note 
that an asymmetry arises in this broadening, particularly 
for low Q, as Eq. El is actually the difference of two con- 
volutions, each weighted by Q'. 

To illustrate this asymmetry, consider the effect that 
a convolution resulting from finite particle sizes has on a 
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FIG. 6: A structure factor peak from spherical particles has 
a peak maximum and peak average below the ideal peak po- 
sition Qo. The asymmetries of the peak are illustrated in the 
inset for Q d = 10, 20, and 40. 



single ideal Bragg peak from an infinite crystal S C (Q) = 
S(Q — Qo), with the effects on all other ideal Bragg peaks 
being accounted for in a piecewise manner. A distribu- 
tion of particle sizes with average diameter d and expo- 
nent n = 3 gives a convolution function (shown in Fig. 



5) 



/de(Q, d) = 2d 



(Qd) 2 



[16 +{Qdf 



(21) 



The convolution of the two [Eq. El produces a broadened 
peak for S(Q), shown in the inset of Fig. HOfor Q d = 10, 
20, and 40. A smaller average particle diameter d (or a 
lower value of the peak position Qo) transforms the ideal 
peak of S C (Q) more asymmetrically than a larger average 
particle diameter (or a higher peak position). The peak 
symmetry is restored for very large values of Qod, as the 
subtracted term in Eq. 1171 becomes negligible. 

Both the peak maximum Q max and the average peak 
position (the normalized first moment) Q avo are also 
shown in the figure as a function of Qod. For most values 
of Qod, Qave < Qmax < Qo, which suggests that not only 
does the peak position shift to lower values of Q, but it 
also diminishes slower to the left of the peak maximum 
than to the right. For the extreme case when Q d < 5.75, 
Qmax quickly goes to zero, and as a result Q avc actually 
increases (it continues to broaden as Qod decreases, but 
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only for positive Q) . This occurs when the characteristic 
width of /de(Qj d) is approximately Qq. 

VI. CONCLUSION 

The predictive power of having an analytical form of 
the PDF of spherical particles has been clearly demon- 
strated. In this case we were provided with structurally 
refined experimental PDF data from a neutron source 
and a distribution of particle diameters observed from a 
TEM. Only the two together provided a means of test- 
ing the presented theory. With care, one could arrive at 
the same predictions by deconvoluting the experimental 
structure factor data of gold nanoparticles using Eq. 1171 
but the refinement analysis that afforded us a comparison 
between theory and experiment suggested a real-space 
treatment. Nevertheless, we believe that the convolution 
functions derived here can complement the peak shape 
analysis already used in modern Rietveld packages, 23 
which often use Lorentzian and Gaussian functions alone 
to describe the peak broadening due to particle size ef- 
fects. Ideally, the analytical form of the peak shapes, as 
derived from the convolution functions, should allow one 
to predict an entire particle size distribution, and not 
just an average particle size, by considering together the 
peak maxima, peak offsets, and peak asymmetries from 
an experimental structure factor. 

It was mentioned earlier that the analysis of spheri- 
cal particles should provide immediate insight into the 
effects of nanoscale domains embedded within a host lat- 
tice. When two structures coexist as uncorrelated do- 
mains, for example, a chemically disordered fee structure 
and a second chemically ordered structure with a small 
tetragonal distortion, the total PDF of the solid can be 
taken as the sum of the individual PDFs, each calculated 



with a unique envelope function suitable to their differing 
domain sizes. The broadening of the Bragg peaks result- 
ing from their finite size, combined with the peak split- 
ting of the second tetragonal phase, may dictate whether 
the chemically ordered domains are below the diffraction 
limit, thus requiring a real-space probe such as XAFS to 
accurately measure their local structure. This analysis 
may be useful in explaining the presence of magnetism 
in NiMn alloys when there is no signature of £lo ordering 
of the material apparent in the diffraction data— 

Finally, the real-space treatment of PDF analysis, and 
the subsequent conclusions made about the structure fac- 
tor, is an encouraging approach to solving problems that 
are often only considered in Q space. For example, the 
idea that a particle's surface might have a different struc- 
ture than its core (internal strain), can be realized in real 
space by considering both an envelope function (for parti- 
cle size) , and a convolution function (for varying strain) . 
Can the structure factor be derived from the two taken 
together, just as it was for the envelope function alone, or 
does this problem require the Q-space treatments already 
considered?— 
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